Phase diagram of the spin-Peierls chain with local coupling 
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We explore the ground-state phase diagram of a Heisenberg spin chain coupled locally to optical 
phonons (bond coupling), using large-scale density matrix renormalization group calculations and 
an extended perturbative analysis. For the quantum phase transition from the spin liquid to the 
dimerized phase, we find deviations from previous quantum Monte Carlo and flow equation results. 
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The interaction of electronic and lattice degrees of free- 
dom in combination with reduced dimensionality can lead 
to a variety of interesting effects, one of which is the in- 
stability of a one-dimensional metal towards lattice dis- 
tortion and the opening of a gap at the Fermi surface 
that was first described by Peierls 1 . A similar effect is 
observed in quantum spin chains, where the coupling to 
the lattice can cause a transition from a spin liquid with 
gap-less excitations to a dimerized phase with an exci- 
tation gap. Experimentally such behavior was first ob- 
served in the 1970s for organic compounds of the TTF 
and TCNQ family 2 , but the topic regained attention af- 
ter the discovery of the first inorganic spin-Peierls com- 
pound CuGe03 in 1993 by Hase et alA In this mate- 
rial Cu 2+ ions form well separated spin- 1/2 chains with 
an exchange interaction that couples to high-frequency 
optical phonons (u> w J), and the phonon dynamics at 
the phase transition is governed by a central peak rather 
than a soft-mode behaviorA^ These features distinguish 
CuGeC>3 from other spin-Peierls systems and sparked the 
interest in a non-adiabatic modelling. 

A good starting point is the study of simplified micro- 
scopic models, which can be build from three ingredients, 
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H — H s + H„ 



(1) 



Here H s — X^-^V^i+i an ^ Hp = ^YliK^i describe a 
Heisenberg spin- 1/2 chain and a set of harmonic (Ein- 
stein) oscillators which are coupled by an interaction 
term H sp . For this interaction we can consider two simple 



forms, 



Hf; B = gu>J2<p\ + bM-S t+1 - S^-Si) , (2) 
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The first type of spin-phonon interaction, , 
been studied with a number of methods, including per- 
turbation theory&i, flow equations^, exact diagonaliza- 
tion£ and DMRGiS. The latter approach identified the 
ground-state phase diagram, but also analytically the 
quantum phase transition from the gap- less to the dimer- 
ized phase is rather well understood: For finite phonon 
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FIG. 1: (Color online) Ground-state phase diagram of the 
spin-Peierls chain with local coupling: QMC and flow equa- 
tion (FEq) results of Raas et al.— ^ compared to DMRG 
data and results of 12th order Schrieffer- Wolff transformation 
(SWT). 



frequency u> the spin-phonon coupling g leads to effec- 
tive spin interactions beyond nearest-neighbor exchange, 
i.e., the low energy physics is governed by a frustrated 
Heisenberg model. As we know from the spin model 



H — ^~](Sj-Sj+i + aSi-Si+2) ■ 



(4) 



frustration can lead to dimerization if the parame- 
ter a exceeds a certain critical value (a c — 0.241167 
in this caseiisiSii 3 .) , and similarly we obtain a finite 

For the second type of spin-phonon coupling, Hl° c , 
which applies to CuGeOs, to date the precise location 
of the phase boundary was arguable. In previous stud- 
ie a 7 i 16 we calculated g c using perturbation theory, a vari- 
ational ansatz, as well as exact diagonalization of small 
systems. This was challenged by results of the flow equa- 
tion method and, in a limited parameter range, by quan- 
tum Monte Carloii^In this article we present unbiased re- 
sults of large-scale density matrix renormalization group 
(DMRG) calculations, and extend our perturbation the- 
ory by several orders in g, high enough to ensure conver- 
gence. The new results are summarized in Figure ^and 
compared to other approaches. More details now follow. 
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FIG. 2: (Color online) Crossing of the singlet and triplet gaps 
for large (right) and small (left) phonon frequency. Note the 
different finite-size scalings. 



Numerical results: DMRG calculations are certainly 
the most precise numerical tool for studying the low en- 
ergy properties of one-dimensional models, but as yet 
have only been performed for the spin-Peierls model 
with difference coupling 1 - Hfp , models with acoustic 
phononsAi, or with XY-type spin interaction^. Wc 
therefore implemented the local spin-phonon interac- 
tion Hl° c , using a high-performance, parallel version of 
the usual two-block finite-lattice algorithm with up to 
1000 states per block. To detect the quantum phase tran- 
sition from the gap-less to the dimerized phase we use the 
established criterion of the level-crossing between the first 
singlet and the first triplet excitation, which was derived 
for the frustrated spin chain 1 -, Eq. (0J, and has success- 
fully been appliecUfl^ in the case of H^ s . For finite sys- 
tems in the gap-less phase the lowest singlet excitation is 
above the lowest triplet, both becoming degenerate with 
the singlet ground state for system size N — > oo. In the 
gapped phase, for N — > oo the lowest singlet becomes 
degenerate with the ground state to form the symmetry- 
broken dimerized state, whereas the lowest triplet main- 
tains a finite gap. Consequently, the two excitations will 
cross at the critical point. Note however, that for small 
phonon frequency uj the relevant singlet excitation can be 
confused with a copy of the ground state plus an excited 
phonon. This frequency range therefore requires rather 
large N for the correct crossing to be detectable in the 
low energy spectrum of the spin-Peierls models. 

In Figure |21 we show the level difference as a function 
of the spin-phonon coupling for various system sizes and 
two typical phonon frequencies. For u> > 0.5 there is 
almost no finite-size dependence of the critical coupling 
<? c , whereas for smaller frequencies the data scales no- 
ticeably. It is tempting to attribute this different be- 
havior to the cross-over from the anti-adiabatic to the 
adiabatic regime, naively expected for u> ~ J = 1. How- 
ever, as was pointed out by Citro et al^ using bosoniza- 
tion techniques, the relevant scale for the adiabatic to 
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FIG. 3: (Color online) Convergence of the Schrieffer- Wolff 
approach with increasing expansion order 0(k). Main panel: 
phase diagram showing DMRG and expansion data for or- 
der 4-12 without and with oscillator shifts. Inset: a e s as a 
function of g for u = 50 and orders 4-12. 



anti-adiabatic cross-over is given by the excitation gap, 
u> ~ A. The larger iV-dependence of g c observed in Fig- 
urc|3for uj = 0.1 is therefore related to the finite-size gap 
still being close u>, and will disappear for N — > oo. Note 
also, that the bosonization results support the analytical 
approach we present in the following, which is based on 
the assumption of anti-adiabaticity. 

Analytical results: Already in our earlier work we sug- 
gested construction of an effective spin model for the 
low energy physics of the spin-Peierls problem by inte- 
grating out the phonons with a Schrieffer- Wolff trans- 
formation^ (which removes direct spin-phonon interac- 
tions) followed by an average over the phonon vacuum. 
In more detail, we apply the unitary transformation 
H = cxp(S)H exp(— S), where 

S = gJ2( b l- h ^-S l+ i. (5) 

i 

Unfortunately, this transformation cannot be evaluated 
exactly, but needs to be approximated by an expansion 
in iterated commutators: 

H = Y i [S,B\h/k\, (6) 

k 

where [S,H] k+ i = [S,[S,H] k ] and [S,H] = H. For 
increasing expansion order these commutators quickly 
become very complicated and easily involve millions of 
terms. Using FORM22S, an algebra tool popular in high 
energy physics, we are now able to push the limit of the 
expansion to order k = 12, a tremendous advantage over 
our previous results with k = 4. For the last step, we 
decided to be more general by averaging the transformed 
Hamiltonian over coherent states with = in- 
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TABLE I: Expansion of the effective long-ranged exchange J„Si-Si+ n in powers of the spin-phonon coupling g. 
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stead of just the phonon vacuum: 

h cS = m\c 

7 

= JoN + J n SfSi+ n + multi-spin terms 

i n—1 

(7) 

This allows for a direct comparison with the flow equation 
result of Raas et alJ£, which is equivalent to a 2nd order 
expansion (k = 2) with phonon shift £ = —(Si-Si+i)g « 
(1/4 — In 2) g w 0.44 g. Here the spin correlator is approxi- 
mated by its value in the isotropic Heisenberg chain. For 
a comparison see the small squares and the thin dashed 
line in Figure ^ 

In Table U we list the expansion coefficients of the re- 
sulting long-ranged exchange interactions J n that con- 
tribute to our effective Hamiltonian Eq. (JJJ. We neglect 
interactions that involve more than two spin operators. 
The phase transition line in Figure^is obtained by equat- 
ing the effective frustration a c s := J2/J1 with the criti- 
cal value a c of the next-nearest-neighbor spin chain (see 
Eq. where for the phonon shift we use the "neutral" 
value £ = 0, i.e., the vacuum. Except for a constant off- 
set, the analytical result matches the DMRG data quite 
well. In particular, the slope of the critical line is cap- 
tured correctly. 

The convergence of our analytical approach is illus- 
trated in Figure |3 depicting the phase transition lines 
obtained for increasing expansion order. We observe an 
oscillatory behavior, such that reliable data results only 
beyond order k > 8. For lower orders, cv G ff may not even 



reach the critical value, as is evident from the order 6 
data shown in the inset. Nevertheless, for order 6 the 
phase transition can be estimated using a Shanks trans- 
formation^. 

We also checked, if the phonon shift £ w 0.44g pro- 
posed by Raas et alJ^ leads to an improved effective 
description, but clearly the dot-dashed line in Figure |3| 
deviates qualitatively from the DMRG data. Within the 
flow equation approach this shift was motivated by the 
technical requirement for normal ordering such that the 
phonons couple to {Si ■ jSj+i — (Si ■ iSj+i)), but not by, 
e.g., some variational principle. Hence, the discrepancy 
is not too surprising. Interestingly, for u> < 3 we obtain 
an almost perfect description of the phase transition, if 
instead we assume a phonon shift of £ « OAAg/uj. As 
yet we did not find a reasonable physical motivation for 
this w-dependence, and the good agreement might be ac- 
cidental. 

To summarize, using DMRG we obtained the, to date, 
most precise numerical result for the location of the quan- 
tum phase transition from the spin liquid to the dimer- 
ized phase in the one-dimensional Heisenberg model with 
local coupling to optical phonons. In addition, we proved 
the convergence of the unitary transformation approach 
that maps the full spin-phonon model to an effective frus- 
trated spin model and allows an analytical calculation of 
the phase boundary in good agreement with the numer- 
ical data. 
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